Non-existence of strong coupling two-channel Kondo fixed point for microscopic 

models of tunneling centers. 



(N 

O 
O 

>.: 



(N 

C/3 



cd 



I 

o 
o 



> 

00 

o 



cd 



o 
o 



X 



I. L. Aleiner ("''') and D. Controzzi^'^' 

'"' Department of Physics and Astronomy, SUNY at Stony Brook, NY 11794, USA 

'''' Physics Department, Lancaster University, LAI 4^3, Lancaster, UK 

^'^' Department of Physics, Princeton University, Princeton NJ 08544, USA 

We consider the problem of an heavy particle in a double well potential (DWP) interacting with 
an electron bath. Under general assumptions, we map the problem to a three-color logarithmic gas 
model, where the size of the core of the charged particles is proportional to the tunneling time, Ttun, 
of the heavy particle between the two wells. For times larger then Ttun this model is equivalent to 
the anisotropic two-channel Kondo (2CK) model in a transverse field. This allows us to establish 
a relationship between the microscopic parameters of DWP and the 2CK problem. We show that 
the strong coupling fixed point of the 2CK model can never be reached for the DWP problem, in 
agreement with the results of Kagan and Prokof'ev, [Sov.Phys. JEPT, 69, 836 (1989)] and Aleiner 
et al, [Phys. Rev. Lett. 86, 2629 (2001)]. 



I. INTRODUCTION 

The possibility to observe the two channel Kondo ef- 
fect in real physical systems has attracted a lot of inteE.- 
est since the pioneering work by Nozieres and BlandinEJ, 
but at present no experimental realization has been con- 
clusively demonstrated. The difficulty lies in the fact 
that the non-Fermi liquid (NFL) fixed point of the two 
channel Kondoj-jmodeH is unstable to various symme- 
try breakinga^'tl which turn out to be important for 
various experimental situations. In particular, channel 
anisotropy is a relevant perturbation and the exact chan- 
nel symmetry is required for the NFL physics to be ob- 
served. In a conventional magnetic realization of the 
Kondo effect this would require the exact degeneracy be- 
tween atomic orbitals that cannot be obtained in any real 
system. In the search for 2CK effect, systems which are 
less directly related to the conventional Kondo physics 
were considered. It was suggested that non magnetic im- 
purity tunneling between two sites and interacting with 
an electron bath could be modeled as a two channel 
Kondo systera. in which the spin plays the role of the 
channel indexQQ (see also Rcf. 10 O for a review). Fol- 
lowing this suggestion, the NFL behaviour of the 2CK 
fixed point was used in Ref. to interpret the low tem- 
perature transport data for narrow metallic constrictions. 
Indeed, in the limit kpa <C 1, where kp is the Fermi wave 
vector and 2a is the distance between the two minima of 
the DWP, only electrons with two spherical harmonics 
[l ~ Q and ? = 1, m — Q) strongly interact with the 
heavy particle. Usually the mapping of the DWP into 
a two level system (TLS), which behaves like a localized 
spin, is achieved by restricting the Hilbert space of the 
atom to the two lowest energy states associated with the 
two minima. (As we will discuss in the following this ap- 
proach is not well justified because it fails to capture the 
connection between the coupling constants of the effec- 
tive theory.) Then the orbital degrees of freedom play the 
role of a pseudo-spin while the real electron spin repre- 



sent the channel index, and, in absence of magnetic field, 
the channel degeneracy was guaranteed by construction. 
A disadvantage of this realization for the 2CK effect is 
that other relevant terms appear which are no longer 
forbidden by symmetry, the most important being the 
spontaneous tunneling of the impurity between the two 
minima. This corresponds to magnetic field in the con- 
ventional Kondo problem and then the resulting Kondo 
model has the form 



H 



Here 



Kondo 
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Hi 
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hS"^ + A,S'. 
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(2) 



represents the propagation of free band electrons, a — 
1,2 is the pseudo-spin index and i = l,2,...,fc is the 
channel index. In the specific case k = 2. The sum over 
repeated indices is implied. The second term in dl^) de- 
scribes the interaction between the localized impurity, S, 
and the electrons 



Hr^l 



Aa^i,,(r)<^#,.(r)5'^ 



(3) 



r=0 



where ct^o are the Pauli matrices in the pseudo-spin 
space, the fermionic operators 4'a,i{f) are the counter- 
parts of ^a,i{k) in the coordinate representation, and v 
is the density of electronic states at the Fermi energy per 
one channel, which is introduced here to make the cou- 
pling constants dimensionless. There are two terms in 
the Hamiltonian (|l|) that break the rotational symmetry 
in the pseudo-spin space. The term proportional to A^ is 
related to the original asymmetry of the DWP, and may 
be greatly enhanced by the disordered potential acting 
on electronsliJ. However, it is still possible to imagine, 
that the original DWP is symmetric, and all the other 
impurities are remote from the DWP, so the value of A^ 



is negotiable and may be even set to zero in some partic- 
ular cases. Nevertheless, the effective magnetic field, h, 
is related to the tunneling of the impurity in the DWP, 
and so is the constant X^ . The relation between them re- 
quires microscopic consideration and one is not allowed 
to neglect h but keep A^, finite. 

The modeL-( l|) is known to scale to the strong coupling 
fixed pointE'E^ where it has a NFL behaviourcl. This 
regime is achieved if both temperature, T, and 



A = v/a1+7i2 



(4) 



do not exceed the Kondo temperature, Tk- Since the 
interaction in highly anisotropic, Aj, can be chosen to be 
equal to zero and A^, <C A^, the Kondo temperature is 
given by (see for instance Ref. fsh 



Tk = D{X.X.f'' ( ^ 



1/2A, 



(5) 



where Z? is a high energy cut-off. In order to check 
whether the 2CK effect can be observed in a DWP sys- 
tem one then has to compare T}~ with. A... Tlusjiroblem 
has recently attracted a lot of attentipnErElllj'EjEj. Some 
alarming results already appeared irJj and were confirmed 
ir£3. The crucial point consists in determining the vari- 
ous coupling constants in ffl) and also high-energy cut-off 
D in Eq. (|5|) starting from a general microscopic descrip- 
tion of the problem. i^i^ 

As first pointed out by Kagan and Prokof'evQa, the 
TLS is not a good starting point since the high energy 
degrees of freedom of the impurity cannot be taken into 
account and they turn out to be essentiaLto-iestablish the 
correct mapping to the Kondo probleniD'B'ta. In particu- 



lar it was shown by model calculation of Ref. 18 that D 
is of the order of the energy distance to the third excited 
level of the atomic system which is several orders of mag- 
nitude smaller than the Fermi energy. Together with the 
restrictions to the coupling constant it leads to h > Tk 
and thus to the impossibility to reach the strong coupling 
limit. However, a question arises, whether those conclu- 
sions are model-dependent or there is a general principle 
for wide class of microscopic models preventing the exis- 
tence of the strong coupling regime. 

In this paper we use a general non-perturbative ap- 
proach to the problem that takes into account all the 
states of the DWP and that, we believe, provides a con- 
clusive answer to the question weather it is possible to 
observe the 2CK effect in the problem of a tunneling 
impurity. We start from a microscopic description of 
the moving heavy particle interacting with electrons in 
a metal. The potential in which the impurity moves is 
chosen to have the most general form with two miniraa. 
(the situation of a potential that presents three minimaEil 
is not considered here). The interaction of the heavy 
atom with the electrons is chosen to be instantaneous 
density-density interaction, which, we believe, contains 
all the essential physics; effects of retardation due to the 



electron screening or the inter-band excitations are small 
either as inverse band gap or plasma frequency and will 
not be taken into account. 

Using a semi-classical approximation to describe the 
dynamics of the impurity and the condition kpa ^ 1 we 
map the problem to a one-dimensional logarithmically 
interacting gas model (LGM). Since a siniilaLmapping 
can be also obtained for the Kondo modeSa, this al- 
lows us to establish a general relationship between the 
microscopic parameters of the DWP and the coupling 
constants of the Kondo model. Our results are in agree- 
ment with the predictions of Rcfs. |3i§7|l^ and show that 
the mapping of the DWP to the 2CK model can be done 
only for energies below ?iT(~„, where Tum is the tunnel- 
ing time of the heavy particle between the two wells. It 
is the existence of this small energy scale, together with 
the relationship between h and X^ , that makes the Kondo 
temperature too small for the strong coupling fixed point 
to be achieved in such systems; Tk turns out to be al- 
ways smaller than A, which makes the strong coupling 
limit non-accessible. 

The paper is organized as follow. In the next Section 
we describe the model and we map it to a logarithmic gas 
model in Sec. III. In Section IV we obtain the relation 
between our microscopic model and the 2CK model and, 
in the last Section, we summarize the results and draw 
some general conclusions. The problem of the effect of 
the electron-hole asymmetry on our results is discussed 
in the Appendix. 




-a ^ 

FIG. 1. Schematic representation of a DWP with two 

minima. The dotted lines represent the lower energy levels 

and A 2 is the level anisotropy. 



II. THE MODEL 

The partition function describing an heavy particle 
moving in a DWP and interacting with an electron bath 
can be written in the following way: 

Z = J D[q{T)]Dm^,T)]Dm^,T)]e-^-i''^^^^l (6) 

where Se is the Euclidean action 

Se^ f drC (7) 



with 



kpa <C 1, 



(16) 



C = C, 



C, 



SC. 



(8) 



Through this Section we will consider only T = 0, an 
extension to finite temperature is trivial and will be dis- 



cussed at the end of Sec. ID 



In Eq. (|^) Cat describes the dynamics of the heavy 
particle and is given by 



^-f(Si +>'")' 



(9) 



where q^r) represent the position of the atom and V{q) 
is the confining potential, which we assume to have 
two minima like in Fig. n, and assume that the level 
anisotropy is absent, A^ = 0, or more precisely, 



V{q) = V{-q). 



(10) 



Note that here and in the following the potential enters 
the Lagrangian with a minus sign with respect to the 
usual definition of the Lagrangian since we work in Eu- 
clidean space. 

The bare electronic Lagrangian has the form 



C„i — 



with 



Ha 



d'^x?/)(x, t)9t^(x, 



d^fc ekV^(k)V'(k) 



Hn 



where ipix.) is the electron field 

V^(x) 



(27r)3 



V(k). 



(11) 



(12) 



(13) 



For the moment wc do not include spin, as it is only a 
spectator. It will be reintroduced at the end where it will 
play the role of the channel index. We chose a general 
density-density interaction between the electrons and the 
atom of the form 

5C = 6H = \ I <i?xij{x, t)V'(x, r)5(x - g(r)) (14) 

which, as explained in the introduction, we believe con- 
tains all the essential physics. 

Given the symmetry of the problem it is convenient to 
expand the electron field in partial waves 



Lm l.m 



l.m 



(15) 



where Y;^™ are spherical harmonics and Rk,i = 2kji(kr), 
with ji{x) = y^TT/2xJi+i/2{x) being Bessel functions. 
^From the condition 



where a is the distance between two minima of the po- 
tential, it follows that only the first two harmonics, I = 0, 
and / = 1, TO = 0, strongly interact with the heavy par- 
ticle. Then we can approximate the electron field as 



V^(r) ^ Fo.o (^) / £^J+{p)RMr) 



-"(;) 



dp 
rJ Jq 27r 



i^-{p)Rp.i{r), 



(17) 



where we have introduced ijj+{k) = ipk,o,o and ip-{k) = 
"0^,1,0- The low energy properties are determined by elec- 
tron states close to the Fermi surface, p -^ kp, so that 
we can introduce fields smooth on the scale 1/kp 



ip±{x) 



P^ki 



^e'P-^±(p). 



(18) 



These fields determine the asymptotic behavior of the 
three dimensional field at rkp ^ 1, thus specifying the 
scattering matrix for the states close to the Fermi energy 
shell. In Eq. (^T^ we suppressed all the higher angular 
harmonics not scattered by the heavy particle. Then, us- 
ing the asymptotic behaviour of the Bessel functions for 
rkp ^ 1 and Eq. (Iq), we can rewrite Eq. ^Tj as 



^(r) ^ iro.o (l) [0+(r)e''=-'- - V+(-r)e 
+ Jyi.o (^) [^-(r)e^'=-'- +^_(-r)e 



-ikFrl 



-ikfrl 



47^7-2 



r \r / 



[?A_ (r)e*'=^'' + V--- (-r)e-*'^^''] , (19) 



where we have chosen the direction of the motion of the 
heavy particle restricted along the z axis. This leaves 
us with an effective one dimensional problem with two 
species of electrons. These are left movers on the entire 
axis as in Eq. ( pj| ) or, equivalently, one can introduce left 
and right movers on the half line. We have chosen the for- 
mer approach for technical convenience. With the same 
accuracy we can linearize the spectrum in the vicinity of 
the Fermi surface: e^ = vp(k — kp) and vp\k — kp\ <C -D, 
where D is the energy scale smaller than the Fermi en- 
ergy. 

Before we proceed, let us note that there have been 
some suggestiontO that the electron-hole asymmetry of 
the original electronic band can lead to an enhancement 
of the Kondo temperature. These effects will not be taken 
into account in this Section, and the corresponding dis- 
cussion is relegated to the Appendix. As we will explic- 
itly show there, they can be introduced perturbatively 
and do not change qualitatively our results. 



For the linearized spectrum, the bare part of the 
fermionic Hamiltonian takes the standard form 



Hq = ivF I dx [ip+dx^+ + i^-dxip- 



(20) 



Substitution of Eq. mt) into Eq. (Mj yields 



4tt 



dpi dp2 

2-K 2tt 



(21) 



V'+(pi)-Rfej=.+j9i,ob(T)] + ^i'-{vi)RkF+pis[q{'r)] 
V'+(p2)-RfeF+p2.o[g(T)] + V2,tl)-{p2)RkF+p2Aq{T)] 



Now we neglect pi_2 in comparison with the Fermi mo- 
mentum and use the condition (16) to expand Eq. (pl|) in 
kpq{T) up to the second order. The smallness of this pa- 
rameter is guaranteed by Eq. ( [l^ ) together with the fact 
that the dynamics of the heavy particle is dominated by 
q{T) ~ ±a. We obtain 



SH = SHo 
6Ho 



5H„ 



\kl 



■4^+i'+ 



kla'' 



'0_^_ 



5Hq = 



Xkl 



kp 



'■-q\T)\ 



(V'+V'+ - V'-V'-) 



(22a) 



(22b) 



from which it is clear that mixing of higher harmonics 
would have been of order 0{{kFaY). Here and in the 
following we use the short hand notation 



V- 



^^(+0) -f v(-o) 



when ever there is no spatial integration. The Hamilto- 
nian ( ^2aD describes the electron scattering on a static 
potential created by the heavy particle smeared between 
potential minima. This term does not have any dynamics 
and may be eliminated by a un itary transformation. On 
the other hand, the term ( |22b| ) describes the excitations 
of the electron system by the moving heavy particle, and 
should be treated carefully. 



To get rid of the Hamiltonian (22a), we perform the 
transformation 



Xkl (. kW\ 
tanS+ = — £- I 1 ^ I , tan(5_ 



Xkl 

ttvf 



''''>" (23) 



Sttvf 



in Eqs. (EQ) and (E2|). After this transformation the term 



Hq + SHo acquires the form (|20|), and Eq. (22b) becomes 

■^-V'+] (24) 



vf 4a 



27rA, 



27rA 



a 



2a2 
2a2 









where the dimensionless constants Ax,z,p are 

2fci?a . 

Az = y=- sm(d+ + d_), 

ttVo 

1 f kFa\ . ,^ ^ , cos^ (5+ + cos^ ^_ 

A, = — ^^ sm05+ + 5_ ± , 

27r \ v3 / cos 0+ cos o_ 

, 1 fkFa\ . _ ^ , cos^ (5+ - cos^ (5- ,„^, 

A„ = — -^ sm((5++^_) i ; . (25) 

" 27r V %/3/ CQs5+cos5- ^ ' 

It follows from the condition ( p^ that 

A„Ap<A, <1, (26) 

independently of the value of the coupling constant A. 
Similar conclusions about the values of the coupling con- 
stants was reached in Ref.Q. 

Introducing the pseudo-spin notation 



^=(Jj), t^ = (^+, ^_) 



(27) 



and choosing the Pauli matrices in pseudo-spin space in 
the following representation 



Ty. — \ ^ Q ) , Ty 



i 



-i / ' ^=" ~ V -1 / ' 



1 



(28) 



the expressions ( pO| ) and (|^ can be rewritten in the 
more compact form 



TT X TJ 

° '^ = * / dxtpd^ip + 27rA^Z(T)t/j^^ (29) 



VF 



2 



2nA,,X{T)iP^xl, + 7:ApX{T)iPiP, 



where we introduced the short hand notation 
Z{t) = -rr-] X{t) = J . 



(30) 



One can see from the Hamiltonian ( p9D that vf can be 
removed by the proper rescaling of the time coordinates. 
In what follows, we will set wf = 1- 

To calculate the partition function it is convenient to 
introduce the U{1) (charge) and SU{2) (pseudo-spin) 
currents 



J = ^p^p, J = ip—ip 

that satisfy the commutation relations 
S'(x) 

[J{x), J(0)]=,^i; 



(31) 



(32a) 



[J,{x), Jfe(O)] = ^^JI!^ + te^'''Six)Ji{xy, (32b) 

47r 

[J{x), Jk{y)]=0. 



The pseudo-spin current commutation relations (32b) de- 
fine the SU{2) Kac-Moody algebra at level 1 {SU{2)i). 
In terms of the currents (^l]) the Hamiltonian ( p9| ) takes 
the Sugawara form 



dx 



-r 



27r ^ 



SH = 2ttA^Z{t)J^ + 2nA^X{T)J^ 



(33a) 
TrApX{T)J. (33b) 



Pseudo-spin and charge degrees of freedom are sepa- 
rated and in the bare part, Hq, one recognizes the 
free charge boson model and the SU{2)i Wess-Zumino- 
Novikov-Witten (W ZNW) model. For a future use, we 
generalize Eq. (33b) by introducing one more coupling 
A„ 



5_H_ 

27r 

y(r) 



= A,Z(r) J, + A,X(t) J, + AyY{T)Jy + ^ 
. Ttun dq 



KpX{T)J 



„ , > (33c) 

2a dr ^ ' 

where i is introduced in front of the first time derivative 
because we work with the imaginary time. The parame- 
ter Ttun has the dimensionality of time and is introduced 
to make Ay dimensionless. It is defined rigorously in the 
next section, however, its exact meaning is not impor- 
tant here. In the original model at high energies Ay = 0, 
however, in our analysis, this coupling will be generated 
in the higher order perturbation theory, as shown in the 
next section. 

As first step in the process of mapping our problem 
to the LGM we need to integrate out the fermionic de- 
grees of freedom. In the framework of the formalism that 
we have chosen here, this requires to construct the Ac- 
tion corresponding to the Hamiltonian (p3). A general 
procedure to do it is provided by non-Abelian bosoniza- 
tion and leads to the WZNW Action (see for instanceE3) . 
For simplicity we chose a different approach that uses 
Abelian bosonization. Despite the fact that the SU{2) 
is a non-Abelian group, the use of Abelian bosonization 
in this case is made possible by the fact that the SU{2)i 
WZNW model has central charge equal to one. Intro- 
ducing two bosonic fields for the charge and pseudo-spin 
degrees of freedom in the standard way 



J = —;=dx 



/TT 
1 



(34a) 



, dx4>s; 
An 

J+ = Jx+ iJy = TT*^ 

Ztt 
D 



J— Jx 



iJ. 



y 



i^/8n(ps 



2tt 



(34b) 



we can write the electronic Lagrangian density in the 
bosonic form 



Ci, = C",+ 6H 



where 



^b^ 2^ [^^-"^P^^^P + i^^^p) 



(35) 



(36) 



and SH can be easily obtained substit utin g Eqs. (Q) into 
( |33c) ) . In the bosonization formulas ( ^4b| ) , D is an arbi- 
trary cut-off required to make the correlation functions 
in the bosonic theory finite. It can be fixed by requiring 
that the correlation function of vertex operators in this 
theory is given by 

(exp(iV4^<?!),(a;))exp(iV4^<?!>,(0))) = -— . (37) 

JJx 

After bosonizing we can integrate out the electronic 
degrees of freedom in (o) and rewrite it as 



D[qir)] e-/^^^-Z4g(r)]Z..,[g(r)] (38) 



where Cat is defined in Eq. (0) and the product 
Zp[q{T)]Zxzy[q{T)] is nothing but the electronic determi- 
nant for the given path of heavy particle ^(t). It has the 
following form 



2p = 



exp [—T^ dTApX{T)J 



= exp 




X(ri)X(T2) \ 



r2 






(39a) 



^-^.-s:(;^)'(2')="7 



J = l 



ndr;dr-i?+(r;)i?_(r7) 



exp 



fl J+(r+)J„(T-) exp ("-4^ I dTA,Z(T)J, 

.7 = 1 ^ -^ 



A. J drZir)±(^-±^^ 



\j=i 



y(- 



JJ dT+ dTj R+ (r + ) i?_ {t- 



j=i 



X exp 



a7-<4(^-^) 



772 ^ rn 






-v)^ 



(39b) 



where 



i?±(r)^^X(r)±^y(r), 



functions Z{t), X{t), and F(t) are defined in Eqs. ([30[), (|33cD , and 

Z, = / exp [-27: f dTA,Z{T)jA \ = exp T^ /"dTidr2 



zin)z{T2) 

(ri-r2)2 



(39c) 



In Eqs. (p9|) averaging over the bosonic fields is defined 



as 



)]e/^-^°. 



DU(x.T 



In the next Section we will use a semi-classical ap- 
proximation to map ( pq ) into the partition function of a 
one-dimensional gas of logarithmically interacting parti- 
cles. This will enable us to establish a non-perturbative 
relationship between our original problem and the Kondo 
model. 



III. MAPPING TO A LOGARITHMIC GAS 
MODEL 



The main goal of this section is to map the model of 
the previous section to the logarithmic ^gas- model. A 
similar approach for the TLS was done incJ'EZl. As usual, 
for such kind of mapping, the first step is to identify the 
logarithmically interacting objects. We will show in Sec. 
[II A that, for the present problem, they are instantim- 
antiinstanton configurations of the tunneling problcmEj. 
In order to find the upper cut-off for the logarithmic gas 
problem V which plays the role of the prefactor in the 



Kondo temperature, we calculate in Sec. IIIB the inter- 
action between the instantons. As we will see, it turns 
out that there is an intrinsic short distance cut-off, that 
makes the theory regular, and is given by the size of the 
instanton. It is analogous to the lattice constant for 2D 
melting problem, or to the vortex core size for supcrfluid 
films. Then a rough estimate for V is given by this quan- 
tity. To define the cut-off with the higher accuracy we 
find the coupling constants corresponding to the chosen 
value of V by microscopic calculation of the energy of the 
logarithmic gas model. The relation between the scale V 
and the coupling constants will be shown to have the 
form of 



h{V) = h [1 + Al (InPn™ + 0(1))] ; 

Ay{V) = -A,A, (iliVTtun + 0(1)) 



(40) 



On the other hand, Eq.(pO|), can be interpreted as the 
result of the first iteration of the renormalization group 
equations for the Kondo problem with Aj, ^ A^,, A^: 



dh{V) 
d\nV 



^hKi 



d\nV 



= -A^A^ 



(41) 



It means that the cut-off D could be chosen arbitrary, 
provided that the coupling constants are adjusted ac- 
cordingly, as it is well-known for any logarithmical prob- 
lem. Accordingly, we will chose the cut-off V to have the 
Hamiltonian of the most simple form: 



Ay[V) = 0. 



(42) 



This equation exclude any ambiguity in the definition of 
V and will allow us to determine the cut-off even with 
numerical coefficient for a wide class of the potentials 
V{q). 



A. Dilute instanton gas approximation 

The partition function (|38| ) can be calculated using a 
semi-classical approximation, which amounts to consider- 
ing only small fluctuations around the stationary points 
of the Action, which correspond to the Euler-Lagrange 
(EL) equations for the problem. A rough condition for 
the applicability of the semi-classical approximation is 
Sinst 2> 1, where Smst is the classical action correspond- 
ing to the trajectory connecting two extrema of the po- 
tential V{q). 

To first approximation we can calculate the stationary 
points using the EL equations for the bare problem 



M^ = VU 



(43) 



which corresponds to the equation of motion for a clas- 
sical particle in the potential minus V[q]. This approxi- 
mation can be improved by taking into account the mod- 
ification of the optimal trajectory due to the electrons, 
however, it introduces only parametrically small changes 
of the coupling constants (because the shape of the clas- 
sical solution is a rigid mode) and does not change any 
of our conclusions. Beside the trivial solution qir) = ±a, 
Eq. (p3|) also admits solutions of the form 



q{T) = ±af{T - Ti) 



(44) 



(kink and anti-kink respectively) with f(±oo) — ±1. 
This solutions are called instantons because they pro- 
duce an almost instantaneous blimp in the Lagrangian. 
The action is invariant under translations of the instan- 
ton center r^, which reflects the time translation invari- 
ance of the original Lagrangian (ph. The instanton is 
characterized by its bare action 



dTCatiafir)], 



and by the tunneling time 

Ttun^ f dT[l-fiT) 



(45) 



(46) 



By construction, the integrand is non-zero only within 
the core of the instanton, so rt„„ has the meaning of the 
size of the core. 

In what follows, we will write the explicit results for 
the model potential 



Viq) = 9{a' - qy 



(47) 



even though our considerations by no means are re- 
stricted for such potential. For the potential (|47|), one 
easily finds 



2, ,3 



^in 






UJT 4 

/(r) = tanh --; rt„„ = -, (48) 

Z U) 



where uP' — %go? /M is the frequency of the harmonic 
oscillation in the extrema of the potential. 




q(t) 



FIG. 2. Schematic form of the instanton solution (U4) 
(solid line). The dotted line represents instantaneous spin 
flip in the Kondo problem. 



The partition function must sum over multi-instanton 
solutions satisfying initial and final conditions. Since the 
instanton core is small, it is a good approximation to 
consider multi-instanton solutions (see Fig. ||) as if in- 
stanton and anti-instantons were dilute (dilute instan- 
ton gas approximation), which requires Jt^ — Tj| ^ -&-„ 
(this assumption is completely justified for Sinst ^ lE3). 
Then, in the absence of the interaction of heavy particle 
with electrons, a general trajectory starting and ending 
at point —a can be written as 



9n(T) 



2n 



l + > X-lF/(r-0) 



(49) 



where the positions of the kinks S,2j+i anti-kinks S,2j are 
subject to constraints 



e. < e 



j+i- 



(50) 



The action corresponding to a configuration with n 
instanton-antiinstanton pairs with the exponential accu- 
racy is 



drCn 



9n(T) 



= 2nSi, 



(51) 



i.e. kinks would not interact with each other if the heavy 
particle were isolated from the electron system. 

In the dilute instanton gas approximation, the low- 
temperature partition function of an isolated heavy par- 
ticle can be rewritten also in a form 



Z, 



' 2n 



n=0 



n 



l/T 



d^o = cosh 



T 



(52) 



Hereinafter, integrals over ^j are calculated with the con- 
straint (^2)- The tunneling splitting between two lowest 
levels of the heavy particle is estimated as 



K Y bi^ 



-S,„ 



1 






(53) 



The numerical factor k conies from the instanton deter- 
minant due to the integration around the saddle pointed. 
For potential (|^, one finds k = ^24/7r « 2.76. 

Closing the subsection, we discuss the limits of the 
applicability of the dilute instanton approximation. If 
the potential V{q) does not have singular points the 
main condition of the applicability is Sinst 3> 1. We 
emphasize, that this condition has nothing to do with 
the number of levels localized in each well. In par- 
ticular, numerical solution of the Schrodinger equation 
with the potential (pTj), shows that the third level in 
DWP lies larger than the maximum of the potential al- 
ready at Sinst < 6.06, (at this point hrtun = 10~^), 
whereas the dilute instanton approximation breaks down 
at Sinst — 2.8, (when the second level crosses the maxi- 
mum of the potential) . 



B. Interaction between instantons 

The purpose of this subsection is to estabhsh the form 
of interaction between instantons and anti-instantons due 
to the dynamics of the electron system. To do so we 
substitute Eq.(^) into Eq. (p8|). The functional integra- 
tion is reduced to the integration around small oscilla- 
tions about the optimal trajectory (p9|). Tlie calculation 
of such instanton determinant is standardE3and may be 
performed without taking into account the electrons, pro- 
vided that the condition (E6h holds. As the result we find 



where Zp^^zy [{C}|"i] ^^'e the integrals ( p9| ) calculated 
over g(T) given by Eq. ([49|). These functions depend on 
the positions of the kinks and therefore produce inter- 
actions between them. In what follows, we analyze this 
interaction in details. 

Substituting Eq. ( ^ into Eq. ( ^ ) and neglecting the 
exponentially small overlap of the instanton cores one can 
rewrite Z{t) and X{t) as 



2n 



Z{r) = ^ -l + E(-l)'/(^-^^-) 



j=i 



2n 






(55) 



It is clear that the function X{t) has peaks at the core 
of instantons and vanishes exponentially otherwise. 

Our strategy now is to substitute Eq. (|55| ) into 
Eqs. ( p9[ ) and perform simplifications using the fact that 
the instant on g as is dilut e. This ca n be easily done for 
the terms (b9a|) and (p9q) while Eq. (39b) requires more 



work. For Eq. (39a) we fin 



interaction between the instantons is not important for 
the dilute instanton gas and we can neglect it. Thus, 
for a configuration of n kink - anti-kink pairs, Eq. (39a) 
takes the form 



Zp = exp(2nA2^S'o), 



(57) 



which is just an independent renormalization of the ac- 
tion for each kink. Such renormalization is nothing but 

the polaronic effect for the tunneling. 

Next, we substitute Eq. ( p5| ) into Eq. (39c) and, omit- 
ting an irrelevant constant term, we obtain 



oo 2n „ 

n=0 i=l -^ J 



CLti dT2 



Z{n)Z{T2) 

(n - T2)2 






Ui{t) = / dTidT2ri2{Ti)ri2{Ti) InVWi - T2 - r 
'?2(t) = ^, SS, = Ui{0), 



(58) 



where energy scale V is introduced here to make the 
argument of the logarithmic interaction dimensionless: 
it does not enter into the expression for the total ac- 
tion. In the following we will see that the natural choice 
is 2? = l/rtum where the tunneling time is defined in 
Eq. ([46|), however, we will keep T) as an independent en- 
ergy scale for pedagogical reasons. 

One can see, e.g. from Fig. g, that function 772 (r) 
decays rapidly at r > rt„„. Therefore, the interaction 
between kinks, Ui{t), is logarithmic at t 3> Tt„„ and it 
saturates at t ~ Ttun ■ This is nothing but the manifesta- 
tion of the usual orthogonality catastrophecfl, where the 
high-energy cut-off is determined by the dynamics of the 
heavy particle. 

Substituting ( p8| ) into Eq. (39c) we find that this part 
of the partition function is equivalent to that of the clas- 
sically interacting gas with n positive and n negative par- 
ticles: 



dT,dT2 ^}^'^^^2^ = 2n6So + E Uo{(^ - ^j); 



(ri-r2)2 
SSo = UoiO); Uo{r) = Re 



dTidT2 ?7i(ri)?7i(T2) 

{ti — T2 — T + iO)^ 



(56) 



One can see from Eqs. (|5q) and ( p5D that Uo{t) decays 
rapidly: Uo{t) ~ (T4„„/T)^at r > rt„„. This short range 



Zz — exp(nA_j(55i)exp(— Til) 



(59) 



To begin the manipulations with the contribution 
(p9q), we first study the simplest, tti = 1, term to illus- 
trate the principle, and then switch to the higher order 
terms. We rewrite the prefactor in Eq. (39b) as 



R+{T+)R^{r^) ^ ^ Viir+ - Qrtjr^ - Q ^ ^ t^{t+ ^.W-ir- Cj) 



{t+ -t--] 



i=i 



{r+-r^Y 



«7^J 



(t+ 



(60) 



^" '7+ (7-+ -^jW-i-r- -S,j) 



E' 



(r+ 



2n 

E 



—2 ' '^iW = -TTVliV T (-l)^^Tt„„772(T) 



where in the last line we used once again the fact that the instanton gas is dilute, and functions ??±(t) decay 
exponentially outside the core of the instantons. The first term in the last line of Eq. (pG) keeps times t± close to 



each other, and, as we will see, will produce the renormalization of the action for a single kink. The second term will 
give rise to the interaction between kinks. 



For the exponent in Eq. ( 39b ) one obtains 
1 1 



drZir) 



E(-l)MC^2fe-r+)-C/2(e,-r_)], 



3 = i 



U2{t)= f dTir]2{Ti)\nV 



Tl 



(61) 



where we introduced the arbitrary cut-off D for the same 
purpose as in Eq. (|58|). Similarly to Ui{t), the poten- 
tial U2{t) is logarithmic at t ^ Ttun and it saturates at 

^From Eqs. ( pOJ ) and (|6l|), we can easily rewrite, within 
the dilute instapton approximation, the m = \ contribu- 
tion to Za;zy asEZi 



-■xzy 



= 2n 



A2 
4 



A2 
y 

4 



L/-^ + ly'-^ + 0{KKyK,) 



+ K,; (62) 



Ix = Re / dridT 



VliTl)miT2) 



^2 ' 



Ttu«R.e / dridr; 



(n - T2 + zO)^ 
V2{ri)miT2) 



4Ki 



(ri -T2 + loy 
2 ^ [a,-(-i)^a,] [a, + (-i)^a. 



2n 



exp <^ A,^ (-1)'= [U2{C^k) - t/2(0fc)] 



fc=l 



where we used the short hand notation 

Sij = ?i ~ Sj 



(63) 



for the distance between the instantons, and ^ means 
that terms involving ^u are excluded. Deriving Ki from 



Eq. (|2|) , we used the properties J dr-qi (r) = r^tm , and 
J dT'r]2{T) = 1. The dimensionless coupling constants 
Ax^y entering into Eq. (p2h are defined as 



A., = A, 



A„ = A„ 



AyA, 



dTrj2{T)U2{T)+0{AlAx) 

dr 



-Vi{T)U2{T)+0{AlAy) (64) 



The factor Ki can be also rewritten in terms of the 
partition function of the classical logarithmic gas. In or- 
der to do so, we relate to any jih instanton the additional 
charge 



H = -1, 0, 1. 



(65) 



As we will see immediately, the physical meaning of the 
"neutral" kink, Hj = is the tunneling of the heavy par- 
ticle without excitations of the electron system, where as 
"charged" kinks fij = ±1 describe the electron assisted 
tunneling. Using the notation introduced in ( |65| ) and the 
fact that asymptotic behavior of the potential U2{t) is 
logarithmic, we rewrite Eq. (p2h as: Ki — Km=i, where 
the more general quantity Km is given by 



Kr 



^(te},^"l 



J2 r{{^l,})e^p{-n2 



2n 



^2 = - E (72 (Cy) h-Ai. + (-l)VjA,] 



r({A^.}) = n 1 

3 = 1 \ 




(66) 



The summation over all the configurations of {fJ.j} is 
subject to the condition of charge neutrality 



2n 







(67) 



and to the mth order perturbation theory constraint 

2n 






(68) 



Equations (p3) and (pq) allow for very natural gener- 
alization to the higher order terms. For a configuration 
containing n kinks and n anti-kinks we wish to keep only 
terms of the order of (nA^)™ and neglect the terms which 
scale like (nA^)™ etc. It amounts to the neglecting in a 
product Y[T=i'^('''t^'^^'''7'^ configurations that contain 
more than two kinks coinciding. We employ this approx- 
imation to extend Eq. (62) for general m. To get rid of 



the combinatorial factor 1/(to!)^, we impose additional 
constraints on the integration r- 
The result then acquires the formE3 



Tj+i, andTj- 



< T. 



i+1- 



'rC™) A2m ™ 1 , 

^ fe=0 



(69) 



where the combinatorial factor takes care of the ordering 
of the paired kinks and the factor K^n is given by Eq. ( pq ) 
with the constraints (67) and (pq). 

The total contribution for 2n instantons is obtained by 
summation of all the orders of perturbation theory: 



Z Z ^ x^v 



m=0 



oo m 

EE 

oo 



(Alnl)'^ {Al/4r-' K^^, 



k\ 



e"A^^^(A^/4)'"i^„ 



(70) 






(75) 



We now substitute Eqs. (^, (|9|), and (|5^) into 
Eq. (p4). As result one finds 



oo ^ 

^ — n r . 1 -^ 



■H 



(71) 



=0 {h} 



where the integration is performed with the constraint 
(pO|), r({^j}) are defined in Eq. (|66|), and the summa- 
tion over charge configurations {f^jj is performed with 
the neutrality condition (pq). Level splitting due to the 
spontaneous tunneling renormalized by the interaction 
with electrons (compare it with (53)) is 



sAlUoiO) sAlUiiO) sAlh 



h — h exp 

\ z, z z 

^/,^ + f^W^ , fA|CMO) £A|4 

~ i 2 2 2 



(72) 



with entries defined in Eqs. (g5|), (p6|), ( |5S| ) and (|62|), and 
s = 1, and we put Ay = 0. 

The energy of the classical logarithmic gas H ~ Hi + 
7^2 is found as 



^ = -E^2(ey) 



fijUi + (-l)VjA^ + 



i-iy+^Ai 



(73) 



where we neglected the difference between the potentials 
U2 and Ui from Eqs. (^ and (|8|) for ^^j much larger 
than the size of the core of the instantons. 

So far we have considered only spinless electrons. The 
real spin is trivially included. We notice that there the 
electron spin commutes with the Hamiltonian, and there- 
fore the fermionic determinant for spin 1/2 electrons is 
factorized onto product of two fermionic determinants for 
each spin. It results in the replacement 



P" — > [?■ 1 

■^-^PyXzy ' y-^p^xzyi 



in Eqs. (pq) and (p4). Such replacement is taken into ac- 
count by introducing an additional "spin" for the kink, 
a^ =|, [. By simple repeating of all of the consideration 
of this Subsection, one finds instead of Eq. (|7l|) 



^ — n r . 1 r _ 1 ^ 



-n 



(74) 



= {p,}, {<T,} 

with the Hamiltonian 



and r(|/^,| ) a re defined in Eq. (pq). The expression for 
the action ( |72| ) should be used with s = 2. This consti- 
tute the result of the mapping of our original problem to 
a one-dimensional classical gas of charged particles, inter- 
acting via the potential U2{r) given by Eq. (|6^)e3. The 
fugacity and the interaction potentials have been calcu- 
lated starting from a completely microscopic theory. 

Energies of the logarithmic gas explicitly depend on 
the high energy cutoff V, see Eq. (61), and so do the 
coupling constants (p4) and (^). It is easy to see that 
the dependence of the coupling constants of cut-off in- 
deed has the form of Eq. (EB) (more formal argument 



that Ay is indeed analogous to the corresponding cou- 
pling constant for the Kondo problem will be given in 
the next section after Eq. (pj)). Indeed, with the help of 
Eq. (Ell), we rewrite Eqs. (pj) and ([72[) as 



A„ 



-A,A, [InPrt, 



a-u 



(76a) 



h = h{l + Al [\nVTtun + ah] + A2C/o(0) + A^/,} ; (76b) 



dTldT2 



'7i(n)f?2(T2)ln 



Tl 



ah^ dT2dT2r]2{Ti)r]2{T2)ln 



Ti -r 



Ttt, 



(76c) 
(76d) 



Because J dTr/i{T) = Tt,m, and J dTri2{T) — 1, the pa- 
rameters ah.y, Ix and t/g can take only numerical values 
of the order of unity not depending of the cut-off V. We 
now define V in accordance with the rule (B3) . We obtain 



V 



1 



nu 



(77) 



which together with Eqs. (|76cD, (^J, (|46|), (|5^), (|5§) 
solves the problem of the relation of the high-energy cut- 
off with the form of the instanton solution of the heavy- 
particle dynamics in the double well potential. Parame- 
ter /i, in its turn acquires the form 



h[l + Al [ah - ay) + A2[/o(0) + All, 



(78) 



For illustration purposes we calculate the numerical value 
of the constant for the potential (47). Using explicit form 
of the instanton solution (Eq), one immediately finds: 



C/o(0) 



12 



C(3) 



ay = tth = 7 In 
" 4 



2e^+^ 



1.461, 
« 0.281. 



where C(3) is the Riemann Zeta function, and C 
is the Euler constant. 



(79) 



0.577 
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Closing the section, we write down the result for the 
logarithmic gas model with the cut-off result. For Aj, = 0, 
Eq. ( pq ) for factors r(|/i-,|) are simplified and one ob- 
tains from Eqs. (^J), ( |75|) 



oo 2n 

^=E^^" E n 



with the Hamiltonian 

2n 



A^e-"" 



\f^i\ 



d£,je 



-H 



(80) 



,(-ir+^AO 



(81) 



Again the integration is performed with the constraint 
( pO| ) and the summation over charge configurations {/ij} 
is subject to the neutrality condition (pZ]). The cut-off D 
is the same cut-off as in Eq. (|37|) , or more precisely 



lim 



In Dt + 27r2Re 



fdr.f 
Jo Jo 



dT2{JxiTl)Jx{T2)) 



0. 



Let us now compare the logarithmic gas models for 
tunneling centers (74), (pG) and for the Kondo problem 
(p3|). First of all, direct comparison of Eq. ( |84| ) with 
Eq. {pm shows that constant Aj^ is exactly equivalent 
to the coupling constants of Kondo model indeed and 
Eq. (|l|) follow. Next, we put Aj, == in Eq. (||) and 
compare the result with Eq. (pO|). We immediately find 
that two models become equivalent upon the following 
identification of the parameters 



The results above will enable us in the next Section to 
make connection between the microscopic parameters of 
the model of tunneling centers and the 2CK model. 



IV. RELATION TO THE TWO CHANNEL 

KONDO MODEL AND IMPOSSIBILITY OF THE 

STRONG COUPLING LIMIT 

A similar nmcedure can be done for the 2CK model 
with real spincja. In our notation Eq. (|l|) acquires the 
form of Eq. (p3|) with ( 33c ) of the form 



D^V ^ 



6H = AttX,S,J, + AttXxSxJx + ^i^XySyJy + 2hSx, (82) 

where we included constant Aj, for the sake of generality. 
We extended the results of Refs. |2^j2^ to include the 
presence of a transverse magnetic field by the repeating 
all the steps of the previous section. The only difference 
is that the spin flips are instantaneous (sec Fig. Q), so 
that the high-energy cut-off is determined by the elec- 
tronic scale D. We obtained the partition function for 
the 2CK model (p2|) in the form (0) as 



with 



n=0 



2n 



E 
{P'j},{'^i} 



d£,je 



-Hi. 



(83) 






2h 



ImjI 



(84) 



and 



Hk 



ondo 



Scnaj [^J'j^J.^ -K (-l)Vj2A3 



+ i^iy+^2XHln\D^,,/h\ (85) 



^tur 



h ^^ h K h 



\r ^ A., 



A, 



e'^yV 



Aa; {hrtun) 



(86) 



where the tunneling time is defined in Eg. ( 46 ) , and nu- 
merical constant ay is defined in Eq. 



for the arbi- 
trary DWP and calculated for model (^) in Eq. (^. 

We note that the equivalence between the tunneling 
impurity model and the 2CK model for times larger then 
Tfun, is non-perturbative in the sense that it is estab- 
lished at any order in perturbation theory. We reiterate, 
that the mapping has to be performed with account of 
all of the excited states of the movable atom. 

/^From the relationships ( p6| ) it is easy to show that 
the Kondo temperature is always smaller then h (Az = 
in our model) and then the 2CK regime can never be 
reached. Indeed, inserting Eqs. ( |8q ) into (||), one has 



¥^--(1: 



ihrtun)\ 7-^-^. (87) 



From the condition (161), and (E5h it follows that 



A^<1, 7>1, -r-<l. 



/^From Eq. (^3|), we have usual relation for the tunneling 
splitting of two levels: 

hTtun < 1- 



Together with Eq. (|87|), this implies that 

TK/h < 1 (89) 

and then the strong coupling regime cannot be reached. 
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V. DISCUSSION 

In this paper we considered a general model describing 
a tunneling impurity moving in a double well potential 
and embedded in a metal. The main motivation for this 
work was to provide a conclusive answer to the ques- 
tion weather it is possible to observe the strong coupling 
regime of the two-channel Kondo model in such a sys- 
tem. In order to answer this question one has to find the 
correct relationship between the microscopic parameters 
of the tunneling impurity problem and the coupling con- 
stants of the effective two-channel Kondo model. Previ- 
ous results showed that the two-level system is not a good 
starting point because all the excited states of the impu- 
rity play an essential role. In order to take into account 
all the excited states of the double well potential problem 
we used a different approach. We mapped the tunnel- 
ing impurity model into a one-dimensional logarithmic 
gas model using a semi-classical (dilute instanton) ap- 
proximation to describe the dynamics of the impurity. 
Since the same mapping can be done for the two-channel 
Kondo model, we obtain a general relationship between 
the coupling constants of the two models. This relation- 
ship is obtained taking into account all the excited levels 
of the heavy particle and is valid to any order in per- 



turbation theory. We demonstrated that, in the effective 
two-channel Kondo model (^ , the values of the coupling 
constant A^, and the transverse magnetic field h, are in- 
trinsically related, and can be never considered indepen- 
dently. This fact, together with the existence of an intrin- 
sic high-energy cut-off in the theory, l/rtum conspire in 
such a way that the Kondo temperature is always smaller 
then the effective magnetic field h. Then the strong cou- 
pling fixed point of the two-channel Kondo model can 
never be reached in this system. This results are valid 
for any form of the double well potential and are robust 
against specific properties of the electron system. 
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APPENDIX A: EFFECT OF THE ELECTRON-HOLE ASYMMETRY 

In this appendix we present the general scheme for the calculation of the parameters of the effective model ( |74| ) 
for a general electron spectrum without involving the linearized approximation (Efl) from the very beginning. Our 
motivation for doing so is to provide framework in which further discussion (quite futile to our opinion) of the role of 
electron-hole asymmetry should be performed. 

Integration over the fermionic fields in Eq. (ra) gives the formal result 



D[i^{^,T)]Dm^,T)]e- 



-SE[q,i),lp] ^ ^-S,i[q{T)]~fdTCat{q,'i) 



where the effect of the electrons on the atom is described by 

~Sei[q{T)]=Trln 



-^-^"(P)-^0^ + ^^P)-^(") 



(Al) 



(A2) 



where p is the quasi-momentum of the electrons, ^(p) = diag[^j(p)] is the spectrum of Bloch electrons, and V{r) 
is the potential of the interaction of the atom|»ath the electrons, which can be of more general form than the local 
interaction (p^). Finally, 0(p) is the standard^ non-diagonal component of the coordinate operator in the basis of 
the Bloch functions and it describes the inter-band scattering due to the atomic potential. In Eq. (A2) and thereafter 
all the energies are counted from the Fermi level. 

It will be convenient for us to express all the quantities in terms of the solu tion of the scattering problem on 
immobile atom. In order to do so, we perform the unitary transformation in Eq. (A2) as 



Trln 



= Tr In 



^iq(r)p 



-jq(T)p 



and obtain from Eq. (|^ 



Sei [qir]] = Tr In [ - 9, + ^pq(r) - ^(p) -V\, V = Vli— + f2(p) 



(A3) 
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Our goal now is to expand Eq. (A3) in powers of q, and relate the expansion coefficients to the coupling constants 
of the low energy theory ( |74| ) . To facilitate such an expansion we introduce the Matsubara Green functions and T 
matrix of immobile atom 

1 

Go(ie„,p) 

G(ie„,Pi,P2) 

L ^ ' J pip2 

T(z£„,pi,p2)= V (l - Go{ien)Vy^ (A4) 

-I P1P2 

with e„ = TTT(2n + 1) being the fcrmionic Matsubara frequency, and their retarded and advanced counterparts 



i£n - i{p) ' 






\ 


1 


P1P2 


= G'o(ie„,pi) 


1 — 


_ien - Up) - V_ 


= 


f(i-Go(z£„)1>)"' 








J 


P1P2 



We note the identities 



Gf ^^(e) = Go(£ ± iO); f^^^\e) = f{e ± ^0). 



Go(i£„,p) - Go{is„ + iu!n,p) = iLL!nGo{ien,p)Go{iSn + «^«,p); 



(A5) 



r(ie„,pi,p2) -T{iSm, Pi, P2) 



(27r)3 
G(«e„, ps) - G{iem, Pa) T(ie,„, pa, P2) 



T(ie„,pi,p3) 



(A6) 



where ujm = 2T:Tm is the bosonic Matsubara frequency. Hereafter the momentum integration is performed within 
the first Brillouin zone. 



Using Eq. (A4) we rewrite Eq. (A3) in the form of linked cluster expansion 



Sei[qir)]^J2 



00 nim) 



S' 



siT^ = Tr 



m— 1 



—iqpG 



(A7) 



where we omitted the term independent of q{T) and all of the multiplications should be understood in the matrix 
sense. 

Before we proceed we notice that in the absence of the impurity potential quasi-momentum p is an integral of 
motion. Therefore, the coupling to p of the force with non-zero Matsubara frequency has no effect independently of 
the spectrum ^(p). Using the fact that 



l/T 



drq — 0, 



we find natural result that 



Tr 



qpGc 



(A8) 



(A9) 



Ai) 



Now we perform the actual calculation of the expansion ( A6). According to Eq. (A7), the first order term vanishes, 



0. Using Eqs. ([A5|) and ([A9|) one finds for the second order term 



LOn af3—x,y,z 



(AlO) 



na/3 ^ _r^ /■ A^p"/TrJ Go(ze„,p)f(i£„,p,p)Go(i£n,p)Go(ie„ + ic^n,p) +c^n ^ -cu,, I 

/d'^n (Pk { - 

1¥~\^W^^"^^^^\ ^o*-**^" ^ *^"' P)^(^^» + *^"' P' k)Go(ie„ + iw„, k)Go(ie„, k)r(je„, k, p)Go(«en, p) 



with 
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qK) = / drq(T)e^'^"^ 
Jo 



(All) 



With the help of Eq. ([A^), one rewrites Eq. ( |A10D 

"""^ ^ -^E/ (sI^s (2l^3 (pV -p"fc^)Tr{G*(p)f (^£„ + zc.„,p,k)G'*(k)f(z£„,k,p)}; 

G* (p) = Go(i£„ + icj„, p) - Go(«£n, p) 
Performing the standard trick with the replacement of summation over e„ to the integration, we obtain 



jja/3 



I f de , s f (Pp (Pk / „ 
- tanh ' ' - 



2r j (27r)3 (27r) 



^ri^K-""'') 



xIm<^Tr G+(p)f(e + zK|,p,k)G+(k)f^(e,k,p) - Tr G-(p)f(e + iK|,p,k)G-(k)f^(£, k,p 



(A12) 



where 



G+(p) = Go(e + iK|,p)-G«(e,p) 
G-(p) = Go(e + ik„|,p)-G^(£,p) 



(A13) 



At that point it is important to emphasize that |q(w„)| decays exponentially at large frequencies for the saddle point 
solution (Eil), |q(a;„)| ex e"''^"''^*"". The fluctuations around the sa ddle p oint are also suppressed at large frequencies 
due to the kinetic energy in Eq. (g). The Green functions in Eq. ( A13|) a re analytic functions of energy except the 
branch cut at the real axis. It allows one to expand over |w„| in Eq. ( [AlSJ ): 

ImG«(£, p)] r^(£, p, k) [lmG^(e, k)] f^{e, k, p) ~ ^ImTr 
ImG^(£,k)lf^(e,k,p 



J__9 

|w„| de 

d 



Tr 



aGo«(£,p)~^ 



T«(£,p,k) 



de 



ReTr 



2 de V '^' ' Q^ 

luiG^ie, p)l d,f^{e, p, k) + [Re9eG^(£, p)! f^'ie, p, k) 



^^"'(^'^)f«(£,k,p) 



Substituting Eq. (A14) into Eq. (AlC) we find 



^ef^^E E q.{^n)qU^n)l\u;n\R"''+U;lQ"'^ 
u)n a^—x.y,z \ 



(A14) 



(A15a) 



Here 



i?"^ = 2^|d£^tanh^|^^(pV-p"fc'')Tr%-e(p)]f«(£,P,k)^[£-e(k)]r^(£,k,p). (A15b) 



and the first term in ( Al5b| ) is precisely the contribution which describes the effect of the gapless excitations - 
orthogonality catastrophe. It is present for the constant density of states. For the spherically symmetric case, one 
can easily recover frome Eq. (Al5b) the exponential factor in Eq. (39c). 

The second term in the expression for the Action characterizes the electron-hole asymmetry. It has an explicit 
expression 

d'^p d^k 



^"-l4^^{'''-"-''''%'^'- 



Q2 ; 



(Al5c) 



Qi = Im 



de 



tanh — Tr 

TT 2T 



e ^ 9Go^(£,p)^^, ^ ^.dG^{e,k)^j, 



de 



-r^(£,p,k)- 



de 



T'He,k,p) 



Q2 = 2^Re / de-^ tanh ^Tri S[e - |(k)]f^(£, k, p) 



5,f«(£,p,k)<5[£-e(p)] 



Re5,G«(£,p) 



f«(£,p,k) 
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where the term Q i depends on the spectrum only and Q2 is due to the frequency dependence of the kinetic coefficients. 
Equation ( A15q ) vanishes for the constant density of states approximation, but it is not so for the arbitrary band 
structure. However, we saw already that it generates the contribution proportional to the higher power of tUn- Thus, 
this term produces non-singular correction to the leading term. The characteristic value of this correction may be 
estimated as ~ l/(Tt„„e*), where e* is the energy scale governing the electron-hole asymmetry. In principle, it migh t 
be estimated from the thermopower measurements or from the first principle band structure calculation of Eq. ( Al5q ). 
One can proceed with the similar expansion in the next orders of perturbation theory to recover the constant 
Kx- The structure remains the same, the first non-analytic contribution in u;„ comes from terms which are kept 
in the constant density of states approximation, and the next contribution has higher power of u;„, and thus gives 
the correction small as l/(Ttu„e*). Thus, the electron-hole asymmetry (if treated systematically) can not produce 
anything except parametrically small corrections to the coupling constants in contradiction to claims of Rcf. M. 
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